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The use of tank-treading as a means of propulsion for microswimmers in viscous shear 
flows is taken into exam. We discuss the possibility that a vesicle be able to control 
the drift in an external shear flow, by varying locally the bending rigidity of its own 
membrane. By analytical calculation in the quasi-spherical limit, the stationary shape 
and the orientation of the tank-treading vesicle in the external flow, are determined, 
working to lowest order in the membrane inhomogencity. The membrane inhomogencity 
acts in the shape evolution equation as an additional force term, that can be used to 
balance the effect of the hydrodynamic stresses, thus allowing the vesicle to assume shapes 
and orientations that would otherwise be forbidden. The vesicle shapes and orientations 
required for migration transverse to the flow, together with the bending rigidity profiles 
that would lead to such shapes and orientations, are determined. A simple model is 
presented, in which a vesicle is able to migrate up or down the gradient of a concentration 
field, by stiffening or softening of its membrane, in response to the variations in the 
concentration level experienced during tank-treading. 



1. Introduction 

Microorganisms such as bacteria and protozoa are able to swim in a viscosity dom- 
inated environment through a variety of strategies. Some of them, su ch as amoebae 

others utilize 



and s ome bacteria, exploit deformations in their mai n body ([B erg 



cilia (jBlake Ill97lt iBlake fc Sleighlll974h or flagella (|Blum fc Hines Ill979t iBerg Il2004h 



976) 



still others, such as cyanobacteria, are able to generate travelling waves on their surface 
(jEhlers et Al. Ill996f h n all cases, contrary to what happens at macroscopic scales, fluid 



inerti a plays no role, and microscopic swi mming is essential l y a lo w Reynolds number 
affair (|Lighthill lll957llChildress Ill98ll ) (see lLauga fc Powers I (|2009h for a recent review 



One of the motivations for the interest in swimming at low Reynolds numbers is its rele- 
vance for the future realization of artificial microswimmers, which would have widespread 
applications in medicine and in the industry. Over the years, various propulsion schemes 



gether, or connected t hroug h immaterial links an d springs; see e.g. iPurcell 



Naiafi fe Golestanian 



Stone fc Samuel 1 119961 : 



200 



Avron et Al. \ (120051) ') and continuous (jLighthill 



Purcell ( 


1977); 


(jLighthill 


1952; 



Ishikawa fc Pedlev 1 120081 ) . In all cases, proper design of a mi 



croswimmcr entails a complex optimization problem, which must take into account limi- 



tatio ns, such as those imposed by the scallop theorem (jPurcell lll977HShapere fc Wilczek 
Il989h . 

Recently, progress in mechanical manipulation at the microscale has allowed to realize 
the first examples of artificial microscopic swimmers ([Dreyfus et Al. 1 120051 : lYu et Al 
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2006t iBehkam fc Sitti I l2006t iTierno et Al] 120081 : iLeoni et Al\ l2009h . At the present 
stage, however, most of such artificial swimmers are driven by external fields and the 
problem of an autonomous power source remains under study. Among the solutions that 
have been taken into consideration, various methods of rectification of Brownian motion 
( Lobaskin et Al. I l2008t iGolestanian fc Aidari I l2009h . and mechanical reactions in the 



swimmer body, induced by inhomogeneity in the environment, e.g. a chemical gradient 
(|Golestanian et Al. ||2005t ' iPaxton et Al. ll2006t iPoolev fc Balazs ll20Ci7h . 

Given the fact that a microswimmer typically lives in a non-quiescent environment, a 
possibility that has been taken into consideration, is to exploit the velocity fields already 
present in the fluid as an energy source for propulsion. Such a swimmer would sail through 
the fluid, by a sequence of deformations induced in its body by the hydrodynamic stresses 
in the external flow. A recent example of such "passive" swimming has been illustrated 
(120101) . based on a discr ete swimmer design similar t o the one considered in 



Olla. 



Naiafi fc Golestanian I (|2004h and IGolestanian fc Adiaril (|2008l ) 

It should be mentioned that passive swimming (at mic roscopic scale s) already exists 
in nature. An example is the Fahraeus-Lindqwist effect ( Vand I Il948l ): a red cell in a 
small artery will deform in response to the flow, in such a way to be pushed to the 
vessel center, thus decreasing its fluid-mechanic resistivity. In analogous way, vesicl es are 



able to migrate transverse to a wall bounded shear flow thanks to tank-treading (jOlla 



1997: Su kumaran fc Seifert 1 l200lL lAbkarian et Al. Il2002h. and similar behavior have been 



obser ved in quadratic shear flows as well ( Olla IbOOoT Coupier et Al. II2OO8 ; Danker et Al 
20091) . 

Tank-treading h as already been taken into consideration as a possible microswimmer 



propu lsion system (jPurcell Ill977t iLehsanskv fc Kenneth 1 120081 ) (see also ITierno et Al 



( 20081 ) for a somehow related approach). It is not too much of a surprise, therefore, that 
the optimal strategy for a passive dis crete micros wimmer in a viscous flow turns out to 



be a discrete version of tank-treading ( Olla 2010h . It is natural to ask what would be an 



appropriate design for a continuous counterpart of this device. We shall concentrate our 
analysis on continuous microswimmers whose basic structure is that of a vesicle. 

One of the motivation s for the present study i s that the e fficiency of a discrete sw i mmer , 
of the kind considered in lNaiafi fc Golestanian I (|2004l) and lGolestanian fc Adiari I (|2008f ). 
is rather low. It is in fact cx aSR/R 2 , where a the size of the moving parts, 8R is the stroke 
amplitude, R is the body size; typically: a/R,SR/R -C 1. In the case of a continuous 
swimmer, instead, a/R=l and the efficiency would become cx SR/R. 

The complicated problem lies in the design of an appropriate control system for such a 
device. In the absence of a control system, a simple vesicle, immersed in a linear shear flow, 
will stay naturally in a tank-treading co ndition, provided th e viscosity contrast between 
interior and exterior fluid is not too high (jKraus et Al. 1l99rf ). A sketch of a tank-treading 
vesicle in a linear shear is provided in Fig. [T] Tank-treading will make such a vesicle 
migrate away from a solid wall perpendicular to the shear plane and parallel to the flow, 
but, unfortunately, no other migration behaviors are possible. Migration towards a wall, 
for example, would require an impossible condition, in which, the tank-treading vesicle 
maintains an ellipsoidal shape, with long axis aligned with the contracting direction of 
the strain component of the shear. For similar reasons, no transverse migration would 
be possible in an unbounded shear flow. The missing ingredient is some mechanism to 
generate internal stresses that counteract the straining action of the external flow. 

We want to explore the possibility that a vesicle be able to generate such stresses 
through appropriate modification of its membrane properties, namely, local stiffening or 
softening of its outer surface. (We hypothesize that the energy required for stiffening and 
softening of the membrane be negligible compared to the one that would expended to 
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Figure 1. Tank-treading vesicle in a linear shear: the membrane circulates around the vesi- 
cle interior preserving the shape and orientation of the object. A vesicle with a homogeneous 
membrane will maintan an overall ellipsoidal shape, with long axis at an angle tp < 7r/4 with 
respect to the direction of the flow. Low internal viscosities correspond to ip ~ tt/4; high internal 
viscosities will lead to ip — ► and then to transition to regimes in which the vesicle rotates more 
or less as a rigid object. 



achieve an identical migration velocity without the help of the external flow; in this sense, 
we speak of a passive swimmer. Although natural at the macroscale, this hypothesis may 
require further justification in the case of microscopic objects). Notice that the presence 
of a mere inhomogeneity in the membrane, say an inclusion, would be insufficient to this 
goal. Such an inhomogeneity would be advected by the membrane flow and would be 
unable to lead to a stationary vesicle configuration. The inhomogeneity of the membrane 
must not vary with time in a laboratory reference frame, which implies that the mem- 
brane elements must continuously change their properties, as they tank-tread around the 
vesicle. 

The mechanism through which local modifications of the membrane stiffness lead to the 
genera tion of stresses in the fluid, closely resembles the one responsibl e for the Marangoni 
effect ( Young et Al. 1959t Subramanian fc Balasubramaniam 2001 ) . Marangoni effects 

have been considered indeed as a possible mechanism of self-propulsion f or dr oplets in in- 

homogeneous environments ( Kitahata et Al. 1120021 : iFurtado et Al. 2008). and lHanna fc Vlahovska 
(2010J) have studied the Marangoni stresses, generated by an external shear induced re- 
distribution of a surfactant over a droplet. 

The philosophy in the present paper is different: the membrane stresses do not act di- 
rectly, through internal flows in the vesicle, to generate propulsion; they rather contribute 
to modify the vesicle shape, and propulsion is achieved through interaction between a 
fixed non-spherical shape and the external flow. 

We shall consider the case of an ideal membrane, so that the only material property 
that must be taken into account is a bending rigidity. The specific behavior we shall be 
interested in is the transverse drift in an unbounded linear shear, already considered in 
" (|2010h . We shall determine the bending rigidity profile that would generate such 



Olla 



behavior, and investigate the possibility that the required bending rigidity profiles be 
obtained as a direct response of the membrane to the external environment, without the 
need of an "intelligent" , internal control system. 

An analytical treatment of the problem is p ossible o n ly in the c ase of quasi-s pherical 
vesicles, and, for this purpose, the analysis in ISeifert I (Il999h and lOlla I (l2000h will be 
generalized to the case of an inhomogeneous membrane. As discussed in lFarutin et Al~\ 
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, the perturbative problem is singular and care must be taken to scale appropriately 
the shear strength and the viscosity contrast between inner and outer fluid, in function 
of the dev i ation from spherical shape of the vesicle. For the sake of simplicity, contrary 
to Seifert ] (Il999h . all finite temperature effect, will be disregarded in the analysis. 

This paper is organized as follows. In Scc.[5]thc bending forces exerted on the ambient 
fluid by a quasi-spherical inextensible membrane, will be calculated, generalizing to th e 
case of an inhomogeneous bending rigidity, the analysis in lZhong-can fc Helfrichl (ll989h . 



In Sec. [31 the shape dynamics of a vesicle with an inhomogeneous membrane, in a viscous 
shear flow, will be analyzed. In Sec. HI the possibility of drift of a tank-treading vesicle 
in an unbounded shear flow will be discussed, and the bending rigidity profiles required 
for drift will be determined. In Sec. [5j a simple model of vesicle, with a membrane that 
changes properties in response to the external environment, will be presented, and the 
migration behavior of the vesicle will be discussed. Section [6] is devoted to conclusions. 
Additional technical details will be presented in the Appendices. 



2. The inhomogeneous membrane 

The shape of a vesicle can be described in terms of the position R(si, S2) of the points 
on the membrane, in function of a suitable set of curvilinear coordinates s\^- For a 
quasi-spherical vesicle, it is natural to work in spherical coordinates r, 9, cf), such that 
R = R(0, 4>)& r , and we write 

R(6,<j>)=R [l + R(6,<f>)], (2.1) 

where Rq is the radius of the sphere with volume equal to that of the vesicle. 

We can decompose the scalar field R in spherical harmonics: R{0,4>) = ^i m Rim 
xYi m (6, (f). In analogous way, vector fields, such as e.g. the displacement 5R(6,(j)) of a 
membrane point initially at R(6*, cf>), will be expanded on the vector spherical harmonics 
basis: 

Ysim = Ylm&r, ^Elm = — , ■ = ■ , Yjyiim = e r X Yez to , (2.2) 

y/l(l + l) 

so that SH = X] M Zm dRftimY film- The basis in (|2.2[) can easily be verified to be orthonor- 
mal: (film\fi'l'm f ) = /Y*, m • Y^i'm' sin0d0d</> = 5 ^ 5 W 5 mm , . 

Assuming that the vesicle has volume V — ^ttRq, it is possible to express the vesicle 
area in terms of the components Ri m through the formula: 

S = (4tt + e)R 2 , e = ^2'( l2 + l - 2 )\ R hn\ 2 + 0(R 3 ), (2.3) 



where Y^'im = J2i>2 J2m=-i ( Seifert 1999() . The excess area e < 1 in Eq. (|2.3[) . which 



parameterizes the deviation from spherical shape, will serve as an expansion parameter 
for the theory. Notice that the I = 1 terms in the sum in Eq. (|2.3[) are identically zero, 
which reflects the fact, that the components i?/=i, m , to lowest order in e, correspond to 
a rigid displ acement of the vesicle. 

Following IZhong-can fc Helfrich I (jl989h . the bending energy of an inhomogeneous 



membrane can be expressed as surface integral 

H B = i J k (2H - C) 2 dS, (2.4) 
where H is the mean curvature of the membrane, which can be written in the form 



( Zhong-can fc Helfrich 
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H 



V 2 R, with n is the unit normal and V 2 the 



Laplace-Bcltrami operator on the membrane; C is called the spontaneous curvature, and 
k is the bending rigidity. For the sak e of simp l icity, we shall assume symmetry of the 



membrane, and set C = 0. Following I Jenkins I ([1977D . to enforce inextensibility of the 



membrane, we include a position dependent surface tension in the energy integral 



n = n B + k / Tds 



(2.5) 



where T plays the role of a Lagrange multiplier coupled to the local area element dS. 

Let us suppose the membrane is able to react to the external environment through 
local variations of its bending rigidity: 



n = Kq[1 + k(R, £)], 



(2.6) 



where, as in iGoulian et Al. I (|1993l ). k is assumed small and will serve, together with e, 
as a basis for a perturbation expansion in Eqs. (|2.4[) . For the moment, we assume the 
profile k(R, t) to be assigned, and postpone to Sec.[5]any consideration on the dynamical 
mechanisms determining its form. 

The membrane will act on the fluid with a force density 



(2.7) 



which will be the sum of a bending force f B and a tension force f T . In the case of a 
homogeneous membrane, the bending force would be directed along the normal to the 
membrane. The space dependence of k produces a tangential force component. In fact, 
the variation of bending energy produced by a deformation field 5H.(9<f>) can be written 
in the form 



sn 1 



6R„ 



2JH 2 W t n 



d6d<p. 



(2.8) 



where Jdddcf) = dS is the surface element of the undeformed membrane, and subscripts 
n and t identify normal and tangential vector components. The tangential contribution 
in Eq. (|2.8j) accounts for the variation of bending rigidity at position R, from tangential 
displacement of a membrane element from position R — (5Rt to R. 

To lowest order in e, i5Rt is a combination of vector harmonics ~Y fJi i m with fi = E, M. 
From the relation fYyilm • VYi' m 'dS' = [see Eq. (|2.2p ]. the tangential component of 
Eq. (|2.8|) has only components from SB-Eim- As it will soon become clear [see Eq. (|3.9[) 
below], this implies that flows on the membrane, induced by inhomogeneity of k, are 
necessarily associated with vesicle deformations. 

To explicitly calculate the bending force, we expand the mean curvature H and the Ja- 
co bian in powers of e. The mean curvature of a quasi-spherical membrane was calculated 



Zhong-can fc Helfrich I (jl989l ). and can be rewritten m the form: 



H 



i-[-l + (1 + iv t 2 )i? - R(l + W 2 t )R + 0(£ 3 / 2 )], 



where Vt = Ro^t- I n analogous way, we can write for the Jacobian: 



J = Rl sing [(1 + R) 2 + \ UdeR) 2 + 

L 2 \ sin 



+ 0(e 3 / 2 ) 



(2.9) 



(2.10) 



Substituting into Eq. (|2.4[) . we obtain the expression for the bending energy: T-L B = 
ko J{2 + 2k[1 - V 2 t R] + RV 2 R + \{^ 2 t R) 2 + 0(e 3 / 2 ) + O(ke)} sin6»d6»d0. Exploiting Eqs. 



6 



P. Olla 



IplB]) and ([2"7F]l . using the expression V 2 = <9 2 + cot6><9 e + (sin fl)^ 2 , + 0(e 1/2 ), and 
expanding on the basis (|2.2j) . we obtain the following expression for the bending force 
density, valid to lowest order in e and k: 

H ° Irn 



- 2y/l(l + l)~K lm Y Elm }6(r - Rq). (2.11) 

To this order of accuracy, the bending force acts on the fluid at the spherical surface 
r = R . Taking for hi m , the spectrum produced by a discrete set of inhomogeneitics 
(e.g. inclusions) in the membrane, Eq. (|2.11[) would lead to the zero tem perature expres 



sions f or the interaction forces among such inhomogeneitics, calculated in lGoulian et Al. 



(|1993[ h Notice that inhomogeneity of the membrane produces a tangential component in 
the bending force, that, in order for membrane area to be conserved, must be counter- 
balanced by tension forces. 

The tension force f T is obtained from variation of J TJd9d(f>. We restrict our analysis 
to a situation in which T = 0(e 1 ^ 2 ) : corresponding to a weak shear regime in which 
the hydrodynamic and the bending stresses, that must be balanced by f T , arc of the 
same order [see Eq. (|3.6p below]. Global area changes are quadratic in R, while local 
changes are linear, hence, it is convenient to separate in the tension, global and local 
contributions: T = T9 l ° + T loc . The variation of the anisotropic part is / T loc 6 JdOd<f) = 
/ STL [2(T loc / R )e r - V t T loc + 0{e)] Jd6d(f>. This leads to the contribution to the tension 
force, to lowest order in e: f T ' loc = -(n a / Rf,)[2(T loc / R Q )e r - V t T loc ]S(r -Rq). 

The isotropic contribution to the tension energy is kqT 91o S; its variation is simply 
K T glo dS, which, from Eq. (|2.3j) leads immediately to the result, expanding in vectors 
spherical harmonics: f T ^° = -(k T3 1 ° / i?g)£ Zm (Z 2 + I ~ 2)R lm Y S i m S(r - R ) (|Seifert I 



19991) . Expanding also f T ' ioc i n vector spherical harmonics, and summing to f T ' 9 '°, we 



obtain 

r = r ,„. + r ,,„ = _ ^ + l - 2)T° l °R lm + 2T l C]Y slr 



lm 

^loc 



- ^T(T+T)T l CYvi m }8{r - Ra). (2.12) 

Setting f t B + f t T = 0, we obtain the local tension T/^ c in the absence of external flow. 
Substituting into the normal component fjf, it is possible to see that the contribu- 
tion by inhomogeneity of the membrane, to the total normal force f n = fjf + f„ , is 
= -^Eimt' 2 + I — 2)ki m Ysim (we shall identify contributions by inhomogeneity 
of the membrane, in general, with superscript "m"). As it could be expected from anal- 
ogous behavior in the case of a homogeneous membrane, the I = 1 components of f m , 
associated with rigid displacement of the vesicle, do not contribute to the sum. Notice 
however that the I = 1 components of f B and i T do not balance in general, and a vesicle 
with an inhomogencous, arbitrarily compressible membrane, could propel itself in a qui- 
escent fluid through a mechanism analogous to the Marangoni effect. A combination of 
inhomoegeneous bendin g rigidity and surf a ce ten sion, rather than just an inhomogencous 
surface tension (see e.g. iKitahata et Al. I (|2002[) ). would be responsible in this case for 



propulsion. More precisely, it is fsim that is directly responsible for propulsion (we recall 
that the induced deformation components SRsim, to lowest order in e, describe rigid dis- 
placements), while /Elm can be shown to generate flows inside the vesicle, analogous to 
the convcction-likc rolls that arc present in a droplet experiencing Marangoni propulsion. 
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3. Deformations in an external shear flow 

We want to determine the shape evolution equation for a vesicle with inhomogeneous 
membrane, immersed in the shear flow 

u(x, t) — 01x2X3 (3.1) 

(.ei = r si n 9 cos 4> : x 2 = r sin 6 sine/) and X3 — rcos9). The derivation closely follows 
the one in Seifert (1999). with additional care, in light of the results in lFarutin et Al 



( 2010T ) . given to the singular behaviors taking place in the limit k, e — > 0. Identifying with 



r]i n and rj out the dynamical viscosities of the fluid inside and outside the vesicle, we can 
introduce dimensionless costants, the capillary number Ca and the viscosity contrast A: 

Ca= 7Wa-Rf and A = ^ Lj (32) 
kq n out 

parameterizing the relative importance of hydrodynamic forces to internal membrane 
stresses, and the ratio between internal and external fluid viscosities. 

The vesicle will produce a flow perturbation u to be added to u on the outside of the 
vesicle, and a flow field u inside the vesicle. The boundary condition at the membrane will 
thus read u(R, t) + u(R, t) = u(R, t), and to this we must add the boundary conditions 
u = at r — > 00 and u = at x = 0. Let us indicate with capital letters values of 
the fluid velocity on the membrane. Expanding in the basis of Eq. (|2.2I) . the boundary 
condition on the membrane becomes: 

In creeping flow conditions, the viscous forces by the fluid are balanced by the reaction 
force exerted by the membrane. To lowest order in e, the force balance at the vesicle 
surface, is evaluated at r = Rq, and the boundary condition equation (|3.3I) is enforced 
at r = Rq as well. To this order in e, from membrane inextensibility, such boundary 
conditions are the ones imposed on the fluid by a rigid spherical surface: 

<L(Ro) = U ( :L = -U«L = -uXm, » - S,E; 



V?« v u; — l-dm — ^ film — "/iZmV 
(°) I n \ — tt(°) — t~t(° 



, lm (Ro) = U^ m = U$> m ee u™ m (Ro), /x = M, (3.4) 

where superscripts indicate order in e 1 / 2 . In the following, although k and R (or e 1 / 2 ) 
are not in general quantities of the same order of magnitude, we shall use superscripts 
to indicate simultaneously order in k and e 1 / 2 [for instance, what we have calculated in 
Eq. (f2~TTj) is actually f B -«]. 

We see that U<?i m = U^l m = U^im = 0. (Besides, it is possible to see that absence of 
external torques implies UMim = to all order in R). 

To lowest order in e, the force balance equation will take the form: 

Uim + aW^rn(U (0) ) - ) + A^im (U (1) )}S(r - Rq) = 0, (3.5) 

where g and g indicate components of the adimensionalized hydrodynamic surface force 
density, associated with fluid flow components vanishing, respectively, at r — > oo and 
r —> (notice that g is defined as as the force density exerted by a flow inside the 
membrane, whence the minus sign in front of the contribution by U (0) ; see Appendix A). 
Inspection of Eqs. (|2.11IA 5HA 6|) shows us that, in order for all the terms in Eq. (|3.5p to 
be of the same order of magnitude, the dimensionless parameters Ca and A must satisfy: 

Ca=0(e 1 / 2 ) and A = 0{e~ 1 ^ 2 ), (3.6) 



corresponding to a regime of weak shear and strong viscosity contrast. 
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The velocity field U determines the membran e dynamics. In particular, the membrane 



displacement rate R(9, </>; t) obeys the equation (jSeifert 11199 



R = U r + U • V t i?. 



(3.7) 



From Eq. (|3T4|) [see also Eq. (|A4 
hence, to lowest order in R: R {1) 



, U (0) 
C/i 1 ' - 



is purely due to the vorticity component of u, 



K=r ' V t-R, where u rot = ^a(x 2 x 3 - x 3 x 2 ), 
that is the vorticity part of the shear flow u. The advection term in Eq. (|3.7|) can then 
be eliminated working in a reference frame rotating with the vorticity of the flow. We 
thus have, for the /j, 



S component of U (1) needed in fflj^U' 1 '): 



rrn(i) 
u Slm 



= R 



lm i 



(3.8) 



and we have introduced a superscript r as a reminder that the components are calculated 
in the rotating reference frame. 

It is possible to see that local inexte nsibility leads to the foll owing relation between 
velocity components on the membrane (jSeifert 1 1199911011a Il2000h : 



rrd) 
u Eln 



JT<1) 

U Slm- 



(3.9) 



Working in the rotating reference frame, Eq. f|3.9[) can then be used, together with Eq. 

(|3.8[) . to express the components UEj m in Eq. ()3.5[) in function of Ri m . 

The /j, = E component of Eq. fj3 . 5|) can be used at this point to express the local 
tension Tf°f in function of the tangential components of f B and of the hydrodynamic 
force. Substituting into the fj, = S component of Eq. (|3.5j) and using Eq. (|3.4[) to express 
U in function of U, we remain with a first order differential equation for the deforma- 
tion component Ri m . Using Eqs. (|2. 11112. 12")) . and (|A 5HA 6[) to explicitate the various 
force contributions, we obtain the equation for the deformation dynamics in the rotating 
reference frame: 



ACaA, di? ^ 



dt 



+ B[R\ m ~ CaCj 



lm 



(3.10) 



where t = at, 



A, 



2Z 3 + 3Z 2 



Cr 
lm ~ 



Di 



1(1 
1 



1) 

4l 3 



6f 



Bi = (I 2 
-41-3, 



Roa 
2(l 2 + l 



1(1 
2). 



1) 



-IL 



■l-2)[l(l- 
21 + 



rpglo-i 



Sim 



Elm 



(3.11) 



That all the 



and the global tension T gl ° is determined from the constrain equation 
terms in Eq. (|3.10l) contribute to the same order in e becomes particularly important at 
the cross-over line in the Ca, A. where the tank-treading regime [the stationar y solution 
to Eq. (|3.10l) ] is only marginally stable. As discussed in Farutin et Al. I ( 201fj) . analysis 
of the crossover region for generic Ca would require inclusion of all 0(e) terms in Eq. 

(|3.10|) . The choice Ca = 0(e 1 ^ 2 ) allows us to circumvent such d ifficultie s. 

For ii = 0, these equations correspond to the ones obtained in lSeifertl ( 1999h and Olla 



( 200(1 ). They differ only for the expression of the coefficient of the time derivative, that, 
to the order considered in e, due to the ordering in Eq. ()3.6j) . must contain only terms 
linear in A. 

In order to return to the laboratory frame, it is sufficient to include in the time dcriva- 
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'Rim 1 -, (3-12) 



dt dt 

Tri' 

Hmm, 1 = u umess m = m ± 1, in which case: 

tilm,m-l = tolm-l,m = ~ m + 1)0 + m )' (3.13) 

Equation (|3.10p takes then the form in the laboratory frame: 

\C&Ai -P- + V" fllmm'Rlm' + BiR hn = CaC hn + Dik im , (3.14) 

L at , J 

m 

where now, from Eq. (| A 4|) : 



Chn = 2iy ^(5/ 2 <5| m |l- (3.15) 

From Eqs. (|3.10[) and (|3.14p . we see that inhomogeneity of the membrane acts in the 
dynamics as a forcing, which acts side by side with the effect of the external flow. Choosing 
components k\ m appropriately, a tank-treading vesicle in an external shear flow could be 
stabilized at orientations otherwise impossible to achieve (e.g. an ellipsoidal shape with 
long axis along the contracting, rather than the expanding strain direction). 



4. Drift behaviors 

Through tank-treading, a vesicle will be able to maintain a fixed shape and orientation 
in a stationary external flow. In the absence of inhomogencitics in the membrane, a tank- 
treading vesicle in the shear flow described by Eq. (|3.ip , will maintain an ellipsoidal shape 
with long axi s somewhere betwe en the stretching direction of the strain and the flow 
direction x% ( Kraus et Al. 19961 ). A fixed orientation is the main ingredient allowing 



migration of a tank-treading vesicle across the velocity lines of the shear flow, and a 
vesicle, in the condition described above, would migrate away from a solid plane wall 
perpendicular to the X2 axis. 

Non-homogeneity of the membrane provides an additional mechanism to control the 
shape and orientation of a vesicle in an external flow, and could be used in principle to 
generate drift behaviors that would otherwise be impossible. We shall focus on the prob- 
lem of generating a transverse drift in the flow of Eq. (|3.ip , in the case of an unbounded 
domain. 

In order for such a drift to be present, it is necessary that the velocity perturbation 
u has components fj, = S, E, / = 1, signaling the presence of a net hydrodynamic force 
acting on the vesicle (see Appendix A). To obtain such harmonics, we must include in the 
boundary condition Eq. (|3.3p determining u, t he effect of non-sphericity of the surface 
r = R{6, 4>). The procedure parallels the one in Olla ( 2000l ). To 0{e 1 / 2 ), we can write: 



LP> = -(i? - R )-?-(u + U W )r=R + U (1 \ (4.1) 
or 



from which we get the boundary condition u^l^ = U (1) . 

Passing to vector spherical harmonics, Eq. (|4.ip will take the form: 



U$ m = - E (nlmlRWm'K,^ + Uj» m> (4.2) 

fi'l'm' 
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where U' = ^(u <0) + u m ) r=Ro . Exploiting Eqs. (|A 2tiA 4() and (|3~4]l . we can write: 



Uki m = 0; 



^Elm — 



21 



-3U< 



Slrn 



2U <0) 



= iV57ra<5 ;2 <5| m |i. 



(4.3) 



It is possible to see that the contribution to drift from Uf$ m vanishes identically. In fact, 
in the rotating reference frame, the components U^V m , fj, = S,E are related to Ri m 



through Eqs. 



and (|3.9p . From Eq. (|3.12p . we find in the laboratory frame: 
2 2 



7T(i) _ r/(i) _ 



it 



Jm 



fi 



Imm 



and we see immediately that U^ m = for /i = S, E; I = 1. 

Returning to Eq. (|4.2[) . we see from Eqs. ()4.3[) and (|A 7|) that the only surviving terms 
in the sum are those for /x = // = E, I = m = 1 and I' = 2, mf = ±1. Expanding 
(film\R\n'l'm') = ^ ; „ m „(^Zm|yi» m "|/Lt / rm')i?;" m », the sum to right hand side of Eq. 
(|4.2[) reduces essentially to two terms, involving matrix elements: 



1 



7 



ll|y 30 |21) = 7A /— , (11|Y 3 2|2,-1> = 



1 



Substituting, together with Eq. (j4~3|) . into ([4T2]| and then into Eqs. (| A 711 A 8|) , we finally 
obtain: 



[7 



drift 



aRo 



5 jjdrift 

—Im[R m ], -L— 
2l7r aito 



5 /7 / 2 - 



21tt V4 V 15 



i? 30 + i?e[i? 32 ] 



(4.4) 



A similar coupling between the I = 2 (or I = 3) harmonics in a shear flow and I = 3 (or 
I = 2) ha rmonics in the intern a l prop erties of a droplet immersed in the flow, has been 
shown in Hanna fc Vlahovska I ( 2010f ) to induce transverse migration of the droplet. 



From Y30 oc cos 0(5 cos 2 8 — 3) and Y32 <x E 2l< ^sin 9 cos 9, we see that a tank-treading 
vesicle drifting to positive X2 will need to have a shape, whose section in the shear plane 
X2X3 (for a > 0) is an egg with the tip at X3 > 0. The geometrical mechanism for drift 
along X2 is illustrated in Fig. [2] an d parallels what is obtained in the case of the discrete 



swimmer discussed in lOlla I (|2010l) . Drift to positive x\ will require, on the other hand a 
shape whose section in the xix 2 plane is an ellipse with the long axis at <f> = 7r/4 with 
respect to X\. A discrete version o f a passive swimmer unde rgoing suc h a kind of chiral 



migration has been illus trated in Watari fc Larson ( 2009f ). In both Olla ( 2010h and 
Watari fc Larson ( 2009( ). the drift was generated in an ensemble of connected spheres 



rotating in a shear flow, imposing a configuration that was on the average asymmetric 
in the laboratory reference frame. 

We can imagine at this point a hypothetical microswimmcr, whose structure is that of 
a vesicle, with full control of the mechanical properties of its membrane, and ask what 
modification of n would be required to achieve the drift behaviors described in Eq. (j4.4[) . 

From Eq. <\AA\ , we see that the drift is maximized if all the excess area is stored in the 
deformation components R^_±2- In order for i?2,±i = 0, we need that the forcing from 
the strain components of u in Eq. f|3 . 14[) , be canceled by the contribution by membrane 
inhomogeneity. From Eq. ()3.15[) : 



1 /IOtt 
«2,±i = t\/ -^-^a, 



(4.5) 



The I = 3 components of k are obtained imposing in Eq. (|3.14p the tank-treading condi- 




a be 



Figure 2. Velocity perturbation (u field) around particles in a viscous fluid, (a): a sphere 
suspended in the shear flow u = 0x2^.3 (a > 0; dashed arrows indicate the strain component of 
u); (6): a tank-treading vesicle suspended in the same flow, shaped to drift towards X2 > [see 
Eq. (|4.4I) ]; (c): a sphere pulled with velocity V to the left. In case (a), 1*2(3:3) = — U2( — £3). In 
case (b), Uz(x3) 7^ — U2(— X3), and ^(x^) has an even component which can be shown to have 
the same sign as the corresponding one in the velocity dipole in (c). In both cases (b) and (c), 
the net effect is a hydrodynamic force pushing the particle to the right. 



tion dR 3 ,±2/di=0, together with R 3m = for m ^ ±2. Using Eqs. (j3lT|) and ([3d3]l : 

^3,±3 = ^-lCaAi? 3i± 2, K3,±l = lCaAi? 3i± 2, 

K3,±a = -10(12 + T^°)R 3t±2 . (4.6) 

We see from the third in Eq. (|4.6p that the component £3, ±2 depends on tension T 9l °, 
which remains undetermined. A lower bound for the amplitude |^3,±2 1 2 can be obtained 
requiring stability of the configuration, i.e. B\ > 0, for which it is sufficient that B2 > 0, 
i.e. T 9l ° > -6. From Eq. ([231), we nnd therefore: 

|£ 32 | 2 > (9/5)e, (4.7) 

and the arbitrariness of K32 reflects the independence of the two expansion parameters 
of the theory e 1 / 2 and k. 



5. A vesicle that "swims" in response to concentration gradients 

A microswimmer, such as the one described in the previous section, would probably 
require a sophisicated control system to achieve the bending rigidity modifications de- 
scribed in Eqs. (|4.6ti4.7|) . One may ask whether a simpler design is possible, in which the 
membrane reacts directly to the external environment, without the need of an internal 
control system. 

We are going to describe such a design, in which the vesicle is able to migrate up (or 
down) a concentration gradient in the shear plane, through softening (or stiffening) of 
the membrane, in response to a local property of the fluid, such as e.g. the presence of 
a chemical substance, a temperature inhomogeneity, or a light intensity gradient. A key 
ingredient will appear to be the presence of a delay in the membrane response to the 
external environment. 

Let us assume that the response of a membrane element in our vesicle, to a concentra- 
tion field O, be described by a linear relaxation equation in the form: 

(d t + U-Vt+7)£ = /?e, (5.1) 

where 7 and j3 could be in general isotropic operators. Equation (|5.1|) could be seen as 
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the result, say, of a process of absorption or chemical reaction with the environment (the 
field may describe e.g. an absorption flux from the bulk, j3 would be an absorption 
constant and oc — V^k may account for surface diffusion effects). The flow on the 
tank-treading membrane is accounted for by the advection term U • Vt^. In stationary 
conditions, the time derivative will drop off Eq. (|5.ip . 

Let us assume the presence of a concentration gradient along X2, so that, on the 
membrane: O = 0(R) = ©o + &R(9, 4>) sm9 sin <j>. (We assume that the diffusive current 
responsible for the gradient in is much larger than its advective counterpart ~ aO' R, 
generated by the flow perturbation due to the vesicle). The effect of the constant part 

00 is an isotropic contribution to h that could be reabsorbed in a renormalization of Kq. 
To determine the anisotropic part, we expand Eq. (|5.ip perturbatively in e: 

(U (0) • V t + 7)« (1) = Q'Rop sin 9 sin cf>, 

(U (0) • V t + 7)k (2) = e'itojftRW sin6» sin - U (1) • V t £ (1) , (5.2) 

and so on to higher orders. 

Let us focus on the regime in which the relaxation time scale for the membrane prop- 
erties is much longer than that of the hydrodynamics, that is a . In this regime, a 
membrane will soften (or stiffen) while cruising at x-i > 0, and start stiffening (or soften- 
ing) when crossing to X2 < 0. In the velocity field described in Eq. (|3.1[) , therefore, our 
vesicle will present a softer (stiffer) side to X3 > 0, and we would expect an egg shape 
with tip pointing at 23 > (at X3 < 0). From Eq. (|4.4j) . this would correspond to drift 
to positive (negative) xi- Unfortunately, a linear theory, based only on the first of Eq. 
(15, 2|) . turns out be insufficient to account for this effect. 

Proceeding as before, we expand Eq. (|5.2[) in spherical harmonics. The lowest order 
contribution to advection reads U <0) ■ Vt = u£=R ' ^> usm § Eq. (|3.12|) . the first of Eq. 
(|5.2|) becomes: 

at'y]fllmm>Kim' + H^tm = Pi@'Ro(lm\ sin 9 sin <p) , (5.3) 

m' 

where now 7; and /3; are numbers. We see immediately that is a superposition of 

1 = 1 harmonics, which do not contribute, to lowest order in e, to the vesicle dynamics 
[see Eqs. (|3. 1013, 1113.141 1], Thus, to lowest order in k, the shape of a vesicle in the shear 
flow of Eq. (|3.1[) will be the same as in the case of a homogeneous membrane: an ellipsoid 
with the long axis between the stretching direction of the flow and the £3 axis. 

At this point, two strategies are possible: one is to replace Eq. (|5.1|) by a nonlinear 
model equation; the other is to take into account higher order terms in e in the vesicle 
response to 0. Now, an equation like (|5.1[) describes the response of the membrane to 
the small variations of that occur on the scale of the vesicle. The physical meaning 
of a nonlinear version of such equation remains therefore unclear. On the other hand, 
the higher order terms in the response to 0, could provide qualitative information on 
the behavior of strongly non-spherical vesicles. This suggests us to opt for the second 
strategy, and to focus on the higher order contributions to the membrane response. 

We must consider the secondary deformations induced by k (2) , and by those non-linear 
contributions to the re-dependent part of the force exerted by the membrane on the fluid, 
that were disregarded in Eqs. (|2.11|) and (|2.12|) . Writing in explicit form: 

P" = f m - L + f m ' w , (5.4) 

with f m > L identifying the contribution by k (2> , f m > N accounting for the nonlinear part of 
the force. 

Notice that the force terms in Eq. (|5.4j) are 0{ke 1 / 2 ). To the same order of accuracy, 
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also O(e) terms should be taken into account; one example are the corrections from 
approximating R = Rq in the Dirac deltas entering i B and f T [see Eqs. (|2.11|) and (|2.12j) ]. 
However, from symmetry of the flow, terms that do not involve k are superpositions of 
even Z harmonics, while drift is produced by Z = 3 harmonics [see Eq. (|4.4[) ]. Such O(e) 
contributions to the force can thus be disregarded. 

5.1. Higher order contributions to the bending rigidity 

Let us consider first the contribution to the membrane force from k (2) , i.e. the term f m < L in 
Eq. (|5.4p . First, however, we have to evaluate the components R^l- From (1, ±1| sin 9 sin (f>) 
= i^/27r/3, and using Eq. (|3.13p in Eq. (|5.3[) . wc find, in the limit 71 /a —> 0: 

K im = 4 y gKdndmo, «= — - — • (5-5) 

In order to determine k (2) , we have to solve the second of Eq. (|5.2[) . From Eqs. (|3.8I3.9|) 

and (EH}, we have U' 1 ' • V t = oEiw -mm^lm-m- Rfi^mm • V. From Eq. (013) 

V H'+rj 

and the fact that k (1) has only Z = 1 components, will be a superposition of Z = 2 
components. Substituting into the second of Eq. (|5.2[) , passing to spherical harmonics, 
and using Eq. (|5.5[) . we obtain therefore: 



\\J\ Yl ^2mw(Zm|Y EW -Vy 10 )^], (5.6) 



where 7; = 7; /a. 

We shall need only the Z = 3 components of k {2) . Solution of Eq. (|5.6p . using Eq. 
(|3.13|) and the expressions for the matrix elements provided in Appendix C, gives then 
the result, for 73 — > 0: 

- (2) - 2K ^lLR' 1 n ) + 2^/Ii^, (1, 



30_ 3x/35V^ 20 ^V3^ 22 



~ (2) / 10 - 5(i) -(2) 2k ~(2)__^ / ^ Su) f 7^ 

K31 ~~ V "7" 21 ' Ks2 ~~ 7/7 K 33 — 3V21 

The contribution to the bending force is in the same form as Eq. (|2.11l) : 

fB ,i„,x = _ 2«o ^ _^ p(i + 1)Yg;m _ ^—^^ _ ^ 

To this, we must add a tension force contribution f T > m > L ; whose effect, as in the derivation 
of Eq. (|3T0)l , is to cancel the tangential part of f B < m < L . Using Eq. (|2T2]| , we obtain 



which will lead to a correction term Dikfz in Eqs. (|3 . 10[) and (|3. 141) . 

5.2. Nonlinear corrections in the bending force 

To lowest order in R, the contribution to force by the Z = 1 components in k is identically 
zero. We must consider terms oc kR in the bending force, which requires keeping terms 
oc kR 2 in the bending energy. The lowest order contribution to the bending energy due 
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to the component is obtained from Eqs. (|2.4I2.9|2.10|) 

1 ^ 



H 



B,in,N _ 



RMV*RV> + ±(V 2 i? (1) ' 



(d e R w )(d 9 R^ 



sin 



| sin 0A9Act). 



Using Eq. (|2.7[) . this corresponds to a force density: 

f B ' in ' N = { [2k (1 » (1 + V t 2 ) V 2 ^ 1 ' + (V 2 £ (1 » ) (V t 2 i? (1) ) 

+ ^ (1) Vf«W]er + 2(V t 2 i? (1) )V t £ (1) } ( 5(r - i? ). 



(5.9) 



As with the other force terms, a tension contribution must be added, that cancels the 
tangential component in Eq. (|5.9p and produces a correction to the normal part, with 

fsi'm' N = V l V +1 /j^'" ,iv [sec Eq. (|2.12|) ]. In terms of vector spherical harmonics, the 
resulting total force will read, using Eq. (|5.5[) : 



in,N 
Sim 



8kq /tt 



m| Viol 2m 



12 



-.(Elm\V t Y 10 \2m)}kR-Z5(r - R ), 



(5.10) 



and again, to determine the drift, only the / = 3 force components will be required. 

5.3. Contribution to drift 

The I = 3 components of the total membrane force are obtained putting together Eqs. 
(|5.7I5.8I5.10[) . Substituting into Eq. (|3.14|) and using the expression for the matrix ele- 
ments provided in Appendix C, leads to the following equation for the I = 3 component 
of the secondary vesicle deformation: 



(5.11) 



where the only non-zero entries of the matrix E 3mm i are: 

1628 



Mm 



E = -*°J± 
3V35' 302 3 V35' 



172 



20 



364y — , E 32 2 — E331 — — — — 



(5.12) 



The lowest order components R^ describe the shape of a vesicle with homogeneous 
membrane, in the flow (|3.1|) . The be havior in a viscous shear flow of a vesicle with 
such characteristics i s well understood ( Noguchi fc Gompper |[2007 : Lebedev et Al. |[2007 : 
iFarutin "a7~I l2010h : the features relevant to the present analysis arc summarized in 
Appendix B. We have: 



nd) 
Jx 22 



J_6d 
V6 



D 1 

H 20 - 



IA 

4 A nr 



where 



A rr = 4 



IOtt 
"3T 



(5.13) 



(5.14) 



Tank-treading propulsion in viscous flows 
0.25 p^- — i 1 1 1 1 1 1 — 



15 




Q I I I I I I I I I 

2 4 6 8 10 12 14 16 
Oct Act- 
Figure 3. Profiles of the normalized drift velocity U drift = -U drift / (aR ke 1/2 ), in function of 
CaA cr for different values of A/A cr . In the four cases: (a): X/X cr = 0.25; (b): X/\ cr = 0.5; (c): 
A/A cr = 0.75; (d): A/A cr = 1. The drift vanishes in the limit A/A cr — > 0. 



is the maximum viscosity contrast for which tank-treading is possible. In the small Ca 
regime considered [see Eq. (|3.6p ], for A > A cr , the vesicle will make direct transition 
to a tumbling regime , in which the vesicle rotates in the shear plane as a rigid object 
( Farutin et ^Tll2010h . 



The membrane tension, entering the relaxation coefficient B3 in Eq. (|5. 1 1[) [see the 
second of Eq. (|3.11|) ]. is determined, to lowest order in e and k, by the shape dynamics 
accounted for by Eq. (|5.13j) : 



T° l ° = -6 + ^-^/l - (\/\ cr ) 2 . (5.15) 

To the order considered in Eq. (|5.11|) . the contribution to tension by the secondary 
deformation R {2> is disregarded. 

Using Eqs. (|5.12M5.15| . Eq. (|5.11[) can be solved in terms of the dimensionless param- 
eters k, ACa and A cr (or e). From here, substituting into Eq. (|4.4|) . the drift velocity can 
be expressed in the form 

Notice that the arguments of JJ dn f t , given the scaling in Eq. (|3.6[) . arc not singular in 
the limit e — s- 0. The profile of the normalized drift velocity U d "f t is illustrated in Fig. 
[3] we see that the maximum is attained for A/A cr = 1 and CaA cr = 0, in which case 

U t ift = ~^\^ 1/2aR o - ~0-25ke^aRo. (5.16) 

0/0UU V 7T 

As expected, softening of the membrane in regions of higher 0, (which implies k < 0), 
will lead to an up-gradient drift of the vesicle. 

The increase of \J dr% ^ t as A/A cr — > 1 and CaA cr —> 0, shown in Fig. [3J is associated 
with a complex pattern of vesicle deformations. Passing from A = A cr to A = for CaA cr 
fixed, the long axis of the ellipsoid described by the components R^m will shift from 
alignment with the flow, to an orientation at 7r/4 with respect to it. At the same time, 
the components R^ with m = 0, 2, that are associated with drift, will go to zero in the 
limit. [This is consequence of the fact that the matrix elements Ei mm , as described in 
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I X 2 

Figure 4. Sketch of the vesicle orientation in the X2X3 plane for A/A cr = 1 and different values 
of CaAcr. From left to right: CaA cr = 0,3.2,6.3,9.5, 12.6. In all cases, k < 0, corresponding to 
drift to X2 > 0. The shapes have been drawn starting from Eqs. (|2.1IB 2[) and (|5.11[) . greatly 
exaggerating the amplitude of the components and R^- 

Eq. (|5.12[) , do not mix even and odd m components, and that the rotation term in Eq. 
(|5.1ip vanishes in the limit. In other words, the fore-aft axis of the vesicle, associated 
with R (2) , and the long ellipsoid axis, align at ir/4 with respect to the flow]. 

For the same reason, the fore-aft and the long ellipsoid axis will align for CaA cr — > 0, 
and X/X cr fixed. At A = A cr , this will occur along the flow direction, which maximizes 
jjmigr ^ n R^ has only m = 0, 2 components). The shift in the vesicle orientation, 
occurring in the process, is illustrated in Fig. 2] 



6. Conclusion 

A microswimmcr could in principle exploit velocity gradients in the suspending fluid as 
a moans for propulsion. We have seen that a vesicle could obtain this result by control over 
the bending rigidity of it s membrane. The addi tional stresses induced in the membrane by 
its own inhomogeneities (jGoulian et Al. 1l993l) could be exploited to counteract the effect 



of the external flow and to allow vesicle shapes and orientations that would otherwise be 
forbidden. 

In the case of an unbounded linear shear, a vesicle with such capabilities could migrate 
transverse to the flow, both in the shear plane and perpendicular to it. Propulsion is 
achieved through modulation of tank-treading by inhomogeneity of the membrane prop- 
erties, which results in a constant fore-aft asymmetric shape for the vesicle (an egg-shape 
with tip along the flow, in case of in-plane drift) and a dipole component in velocity 
perturbation around it, associated with the presence of a transverse force. 

These calculations can be easily extended to the case of a shear flow bounded by a wall, 
thus allowing to consider regimes in which a tank-treading vesicle could be made to drift 
towards a wall, rather than away from it, as occurring normally with an homogeneous 



membrane (jAbkarian et Al. 1120021 ) 



We have shown that such migration behaviors (at least, those confined to the shear 
plane) could be achieved simply by a stiffening (or softening) of the membrane, in re- 
sponse to variations of some local fluid properties, e.g. a chemical concentration. With an 
appropriate choice of parameters in the membrane response, a vesicle could be made to 
migrate automatically, up or down a concentration gradient in the shear plane, without 
the need of complicated control systems. 

The results discussed have been obtained through a perturbative calculation in the 
case of a quasi-spherical vesicle. In this regime, the migration velocity is smaller than 
the velocity scale aRo (the velocity difference in the shear flow at the vesicle scale), by 
a factor ke 1 / 2 , with e the normalized excess area and R the degree of inhomogeneity of 
the bending energy [see Eqs. (|2.3|) and (|2.6j) ]. However, as the mechanisms for stress 
generation by membrane inhomogeneity, and for migration by asimmetry of the vesicle 
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shape, are essentially geometric, the same behaviors are expected to hold also for strongly 
non-spherical vesicles and strong inhomogeneity in the bending rigidity. 

In the present analysis, no attention has been given to the problem of the material that 
could be utilized to synthesize a membrane with the properties described in Eqs. (|4.6I4.7[) . 
or in Eq. (|5.1|) . In principle, the required local modifications of the membrane could be 
achieved by some sort of absorption process or chemical proc ess with the bulk, in an alogy 
with what is done in experiments on the Marangoni effect ( Kitahata et Al. I 2002 ) and 



osmophoresis dNardi et Al. Ill999h . It remains to be determined, however, which mate- 
rial could satisfy the requirements imposed by Eq. (|5.1[) , in particular, slow relaxation 



of membrane properties, compared to hydrodynamic time-scales. At the present stage, 
therefore, as in the case of most microswimmer designs, also the present vescicle-based 
one remains at the level of a purely theoretical model. 

It must be said that the adoption of a vesicle based microswimmer design, rather than 
a microcapsule based one, has been motivated purely by simplicity considerations on the 
constitutive law for the membrane. After all, tank-trea ding behaviors, similar to those of 
vesic l es, are observed also in the case of microcapsules (IBarthes-Biesel 1980: iPozrikidis 



l2003t ISkotheim fe Secomb"ll2007l ). This could open additional possibilities as regards the 
issue of practical realization. An interesting question is whether a microswimmer design, 
utilizing an elastic, rather than fluid membrane, would be more effective in converting 
fluid stresses into migration. 



The author wishes to thank Alexander Farutin and Chaoqui Misbah for interesting 
and helpful discussion. 



Appendix A. Lamb representation 

The Stokes equation for an incompressible fluid are 

t?V 2 v = VP, V • v = 0, 



(Al) 



where P is the pressure. This equation can be solved in the basis of Eq. (|2.2[) , with 
boundary conditions v = V imposed on a surface r = Rq\ the resu lt is the so called 
Lamb representation of the Stokes equation ( Happel fc Brenner 1973f h The solutions for 
r < Rq and r > Rq read respectively: 



"Sim 



J Elr, 



Mir, 



(l + 3-(l + l)y 2 )V slm + jl(TTT){y 2 - m 



Elm 



,1-1 



(i + : 

V M lm V 



(1 - y 2 )V slm + (-(l + !) + (/ + 3)2/ 2 )Fi 



Eh, 



J-l 



(A 2) 



and 



J Slm 



J Elr 



J Mlm 



(I + (2- l)y- 2 )V S im + y/l(l + 1)(J - y- 2 )V E im 



(2 - OVfftC 1 - y~ 2 Wsim + (2-l + ly-^Vi 



Elm 



(A3) 



y 



where y = r/R^. 

The shear flow u = aX2'x.3 and u, that is the flow inside the vesicle, are both in the form 
of v m , while the perturbation u is in the form of v out . The lowest order solutions u (0) 
and u (0) are given by Eqs. (|A3|) and (|A2|) with V = U (0) and V = U (0) . To higher order 
in R, different harmonics in U, U and u, u get mixed due to deviation from spherical 
shape. 
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The vector spherical components on the surface r = Ro for the shear flow u are 
obtained from Eq. (|A 2[> : 



rr(0) 
U S2,±1 



rr(0) 
U E2,±1 



7T(0) 



ai?o V 15 ai? V 5 ai?o 

It is easy to see that ±1 is responsible for the vorticity part of the velocity u r 

50(12X3 - x 3 x 2 )- 

The force density on the surface r 
will read, in dimensionless form: 



(A 4) 



Ro, produced by the internal field in Eq. 



ffSim(V) = ^ 
SEim(V) = ^ ^ 



7 vSlm + 



T/ 



(21 + 1)VW 



(A 5) 



and the one from the external field, Eq. ([A 3 

5Eim(V) = ^ 



2i-+3i+4 T/ 
3a/ TTT^Shn 



i+1 



(2/ 



l)^Eim 



<?MZm(V) = — Q^-VMlm- 



(A 6) 



(A 7) 



The total hydrodynamic force F on the vesicle is obtained from the \i = S, E, I = 1 
components of the outer Lamb solution (|A 3[) u. In particular, for the force components 
along Xi i2 : 

Fi = 2V0^VextRo ReiUglj + y/W^), 

where we have used the fact that, from Eqs. (|3~4|) and (|A4|) : = for (j, = S, E. In 
the absence of external forces, this is cancelled by the drag force —DoU drt ' t on a vesicle 
migrating with velocity JJ dr ' l f t j where, to lowest order in e, Dq is the Stokes drag by 
a spherical vesicle. Now, a spherical vesicle with an inextensible membrane will behave 
with respect to the fluid as a rigid object [except for solenoidal flow components Um/to 
on the surface, that do not couple with U^™-' 7 '; see Eqs. (|3.8[) and (|3,9[) ]. Hence, Do is 
the drag coefficient of a rigid sphere. The drift velocity of the vesicle will be therefore, 
to lowest order in e: 



V drift _ y/Dq, 

where Dq = 6irrj ex tRo is the Stokes drag for a rigid sphere of radius R$. 



(A 8) 



Appendix B. Dynamics of the homogeneous membrane in a shear 
flow 



Assuming a tank-treading regime, the shape of a vesicle in the shear flow (|3.1[) is 
obtained from Eq. (|3.14p setting the time derivative equal to zero. Using Eqs. ()3.11I3.13[) 
and (|3.15p . and setting to zero the inhomogeneous contribution Diki m : 

iCaA ~, 



2 



mi 



/ton 



BoRk 



0. 



IOtt 
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iv6CaA 

R { H + B2R20 = 0, A = A 2 \, 

which gives the result 

= KC&Ks/Q, R%1 = 2KB 2 i, = KC&A, 

where 

Ca 



19 



K = 



IOtt 

B\ + (CaA) 2 V~T' 
Substituting Eqs. (|B 2|) and (|B 3|) into the area constrain (|2.3|) . we obtain 



So 



1607rCa^ 
3e 



(CaA) 2 , 



(B2) 
(B3) 

(B4) 



from which we obtain Eq. (|5.15|) . Exploiting the first of Eq. p. lip , we see that tank- 
treading is possible for 



, . 24 /IOtt 
A<Ac ^23V^T' 



that coincides, to leading order in e, with the result in lFarutin et Al. I(|2010l ). Substituting 
Eq. (|B4| into Eqs. ([B3| and (|B2]| . we obtain Eq. (f57T3|> . 



Appendix C. Matrix elements involving scalar and vector spherical 
harmonics 

We provide below expressions for the matrix elements entering Eqs. (|5. 515. 1115. 12p . 



(30 
(31 

(33 
(31 
(30 
(32 



Y21 sin 9 sin <j)) = 
Y22 sin 9 sin <f>) = 

Y22 sin 9 sin cf>) = i\ 

Y E2 i • W10) = -2 



70' 



70 



14' 



35tt' 



F 10 |20) = 

y 10 |22) = 



3 / 3 
2 V 35^ : 

2 V 7tt' 



(E31|W 10 |21) = 2 



35tt' 



(31|i20 sin ^ sin 0) =i 



35' 



(32 1 y 2 1 sin 6» sin 0) = 

V 7 



(30| Y E20 • VY 10 ) = -3 
(32|Y E22 • VF10) = - 
(31|y 10 |21) = 
(E30|Vr 10 |20) 
(E32|Vyi |22) 



3 

357' 
77' 



6 

357' 



35?r 
1 



7tt 
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